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ABSTRACT 

Cosmological simulations have shown that dark matter haloes are connected to 
each other by large-scale filamentary structures. Cold gas flowing within this "cosmic 
web" is believed to be an important source of fuel for star formation at high red- 
shift. However, the presence of such filamentary gas has never been observationally 
confirmed despite the fact that its covering fraction within massive haloes at high 
rcdshift is predicted to be significant (~ 25%), (Dekcl et al. 2009). In this work, we 
investigate in detail whether such cold gas is detectable using low-ionisation metal ab- 
sorption lines, such as CII A1334 as this technique has a proven observational record 
for detecting gaseous structures. Using a large statistical sample of galaxies from the 
Mare Nostrum N-body-f-AMR cosmological simulation, wc find that the typical cover- 
ing fraction of the dense, cold gas in 10^^ Mq haloes at z ~ 2.5 is lower than expected 
(~ 5%). In addition, the absorption signal by the interstellar medium of the galaxy 
itself turns out to be so deep and so broad in velocity space that it completely drowns 
that of the filamentary gas. A detectable signal might be obtained from a cold filament 
exactly aligned with the line of sight, but this configuration is so unlikely that it would 
require surveying an overwhelmingly large number of candidate galaxies to tease it 
out. Finally, the predicted metallicity of the cold gas in filaments is extremely low 
(< 10^"^ Zq). Should this result persist when higher resolution runs are performed, 
it would significantly increase the difficulty of detecting filamentary gas inflows using 
metal lines. However, even if we assume that filaments are enriched to Zq, the ab- 
sorption signal that we compute is still weak. We are therefore led to conclude that 
it is extremely difficult to observationally prove or disprove the presence of cold fila- 
ments as the favorite accretion mode of galaxies using low-ionisation metal absorption 
lines. The Lya emission route looks more promising but due to the resonant nature of 
the line, radiative transfer simulations are required to fully characterize the observed 
signal. 
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1 INTRODUCTION 

How galaxies get their gas is a long-standing issue. For 
decades, the standard theoretical picture of galaxy forma- 
tion has stipulated that all gas accreted into dark mat- 
ter haloes is shock heated be fore i t rad i atively cools and 
settle s into a g alactic disk (ISilkl 1 19771 : iRees fc Ostrikeij 
Il977l . although iBinnevI Il977l first suggested this need 
not be the case). This picture has recently been revis- 
ited both by analytic studies (|Birnboim fc Dekell l2003l : 
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iDekel fc Birnbpim 20061') and hvdrody naniical simulations 
(both in ID: iBirnboim fc Dekell |2003| and in 3D within 
an explicit cosmological context , Keres et al.l |2005| . l2009l : 
lOcvirk. Pichon. fc Tevssieill2008l : hereafter OPT08). These 
studies have established that in haloes below a critical mass 
shocks are unstable and cannot propagate outwards, so that 
cold diffuse gas and/or cold filaments can penetrate deep 
into the halo without experiencing shock-heating. In con- 
trast, at the other end of the mass spectrum, very massive 
haloes easily sustain a virial shock that is stable against gas 
cooling so that diffuse/filamentary gas is shock heated to 
the virial temperature of the halo as it enters. Finally, at in- 
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termediate halo masses, either gravitationally shock-heated 
hot gas and/or a hot galactic wind coexists with cold in- 
flowing filaments (OPT08): some dense cold filaments are 
stable against the pressure force exerted by the hot gaseous 
material. In other words, the vast majority of galaxy-size 
host halos, especially at high redshift {z > 2), are predicted 
to be threaded by cold gas filaments in a ACDM model 
of structure formation. Therefore, the question which natu- 
rally arises is whether or not the existence of these cold gas 
filaments can be observationally confirmed, and by which 
technique. 

Since the advent of high-resolution spectroscopy has 
enabled astronomers to study the kinematics of the inter- 
galactic medium at high redshifts in a wealt h of detail (e.g. 
Pettini et all I2OO2I : lAdelberger et al.l l2003l : IShaplev et aP 
20031 ). and given the fact that cold gas is thought to flow 
into massive haloes (lO^^M©) at z = 2.5 along fil aments 
with velocities of > 200 kms"^ (jPekel et all l2009l ). it is 
sen sible to think that the spectra of Lyman-break galax- 
ies (|Steidel et al.|[l99^ ) might reveal these filaments as red- 
shifted absorption features. Interestingly, a recent study re- 
ported that few Lyman-break galaxies (LBGs) show red- 
shifted metal absorption lines, suggesting that inflowing 
■as in haloes with 4 x lO" < Mvir < IO^^Mq is rare 
Steidel et a"l]|2010l ). Taken at face value, this appears to 
contradict the theoretical prediction that cold fllaments are 
prominent in the vast majority of high-z halos. 

However, there exist a variety of reasons as to why these 
filaments should be very difficult to detect. The first of these 
is the covering frac tion of the filaments. This is estimated in 
jDekel et al. I I2OO9I ') to be around (~ 25%) for four massive 
halos with Mvir ~ 10^^ M© in the Horizon-MareNostrum 
simulation, counting only relatively dense (A^h > lO'^^cm"'^) 
and cold (T < IQpK) gas within a radial distance 20 < 
r < 100 kpc from the central galaxy hosted by these halos. 
Whereas 25% is indeed a non- negligible covering fraction, its 
exact dependence on reds hift and halo mass remain s to be 
determined. For instance, iFaucher-Giguere fc Kerej (|2010l ') 
also measured the covering fraction for a Milky way-type 
progenitor LBG in a smoothed particle hydrodynamic sim- 
ulation but found a significantly smaller value (~ 2%). Sec- 
ondly, low-ionisation lines can be produced not only by cold 
filamentary gas but also by a galaxy's interstellar medium, 
so that when using a single galaxy to probe the circumgalac- 
tic medium, distinguishing absorption produced by filaments 
from that produced by the ISM is key to prove/disprove the 
presence of cold filaments. Thirdly, there is the geometry of 
the accretion: in order to produce a strong absorption sig- 
nal, a filament needs to be well aligned with the line of sight 
to maximise its column density. Finally, there is the issue 
of metallicity: a metal-poor filament will be transparent to 
metal line observations. 

The aim of this letter is to quantify the aforementioned 
effects to better assess the detectability of the absorption 
signal produced by cold filaments. We show that the ac- 
tual probability of detecting such flows with metal lines 
is much sm a ller th an the high covering fraction derived in 
iDekel et al.] (|20Q9l ') would suggest due to i) the low density 
and ii) low metallicity of the filaments compared with the 
densities and metallicites of the interstellar medium of the 
host galaxy. 



2 SIMULATIONS 

To investigate the statistical properties of cold filaments, 
we analyze the Horizon-MareNostr um simulation, car - 
ried out with the WMAPl cosmology (|Spergel et al.ll2003t) 
using the adaptive mesh refinement code RAMSES ( Tevssien 
2002). Details of the simulation can be found in OPT08, 
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briefly describe the simulation setup and modelling strat- 
egy here. 

The simulation follows the evolution of a periodic cos- 
mological volume of comoving side length L — 50/i~'^Mpc, 
and contains 1024'^ dark matter particles with mdm ~ 
1.4 X lO'^M©. Dark mat ter halos are identifi ed using the 
AdaptaHop algor ithm (|Aubert. Pichon. fc"Colombit ,2004 
iTweed et al] |2009| ). resulting in 3419 halos with Mvir > 
1O"M0 at z = 3.8 and 6456 haloes aX z = 1.5. The spatial 
resolution of the simulation is kept fixed at around 1 /i~^kpc 
physical over the entire redshift range. 

Gas in the simulation can radiatively cool by ato mic 
processes down to 10* A" (|Sutherland fc Dopital ligQal ). A 
fraction of the cold and dense gas (nn > O. lcm"'^) turns int o 
stars following the Schmidt-Kennicutt law (|Kennicuttll998l ). 
and massive stars explode as type II supernovae, redistribut- 
ing energy and metals into the interstellar/intergalactic 
medium. The Sedov Blast wave solution is a dopted for su- 
pernova explosions (|Dubois fc TevssieJbOOSl ) , and reioniza- 
tion is implemented by ins tantaneously turning on a uniform 
UV background at z=8.5 (|Haardt fc Madaulll99d ). 



3 RESULTS 

3.1 Covering Fraction of Dense Gas 

OPT08 showed that the transition mass (Mstream) separat- 
ing cold from hot dominated accretion increases with in- 
creasing redshift. According to these authors, hj z ^ 2.5, the 
cold streams feeding massive haloes (M200 si 3 x lO^^Mo) 
begin to disappear, but based on measurements of a hand- 
ful of ha loes in the Horizo n-MareNostrum simulation at 
z ~ 2.5, iDekel et al.1 (|2009t ) find that some massive haloes 
with M200 ~ IO^'^Mq still show high covering fractions 
(< 25%) of dense {Nn > lO^^cm-^) and cold (T < lO'^K) 
infiowing gas within 20 < r < 100 kpc. In order to obtain 
statistically significant results on the covering fraction in the 
Horizon-MareNostrum simulation, we analyse a// the ha- 
los in the simulation and show the results in Fig. [1] Using 
the complete sample, we find that the average covering frac- 
tion of haloes with M200 ~ lO^^M© at z = 2.5 is ~ 5%, 
about a factor 5 le ss than the covering fraction reported in 
iDekel et al.] (|2009l ) for their sub-sample of M200 ~ 10^^ Mq 
halos at z ^ 2.5. We note that thi s lower value is more con- 
sisten t with the recent findings of iFaucher-Giguere fc Kerej 
(|2010D . However, whilst low covering fractions (~ 5%) are 
computed for most halos at redshifts z < 3, those of higher 
redshift (z ~ 3.8) halos are much larger (~ 25%), as can be 
seen in Fig. [T] This strong redshift evolution in the covering 
fraction of cold filaments between z ~ 3.8 and z ~ 2.5 re- 
flects the aforementioned rapid transition from cold to hot 
dominated accretion. We note that these results are con- 
sistent with the OPT08 value of ~ 10"Mq for Afstrcam at 
z = 4. 
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Figure 1. The covering fraction of cold (T < lO^K) and dense 
(iVu > 10^'') gas within 20 < r < 100 kpc physical as a function of 
the virial mass of halos (Af20o)- Different colours indicate the cov- 
ering fractions at different redshifts. Solid lines include the con- 
tribution from the interstellar medium of satellite galaxies to the 
covering fraction. To exclude this latter contribution, we also plot 
the covering fraction with upper density cut (10^° < A^h < 10^^ 
cm~^, dotted lines). Error bars correspond to the interquartile 
range (25 < / < 75%). For a given halo mass, halos at higher 
redshift show larger covering fractions. This can be understood 
in terms of the rapid development of a virializcd hot medium 
between z = 3.8 and z = 2.5. 

An interesting feature present in Fig. [1] is that the cov- 
ering fraction is higher in more massive haloes at a given 
redshift. This seems to contradict previous findings (OPT08) 
that it is the small haloes which are mainly fed by cold mode 
accretion. It should be noted, however, that the covering 
fraction shown in Fig. [T] does not account for the accretion 
of cold, diffuse gas (cold gas with lower column densities) 
which is only present in halos with masses incapable of sus- 
taining a virial shock at all. Indeed, these small haloes are 
usually located within (or around) filaments whose density 
is low, whereas the filaments around more massive haloes 
tend to be denser. Furthermore, in massive haloes, satellite 
galaxies contribute more importantly to the covering frac- 
tion, as do extended, warped, galactic disks and dense gas 
bridges which result from tidal interactions between galax- 
ies. These latter effects partly explain the trend, but the 
primary driver of the covering fraction increase with halo 
mass is the density of the accreted gas, which is higher in 
more massive haloes. Fig. [T] (dotted lines) substantiates this 
claim by showing how the covering fractions drop when the 
very dense gas (Nh > lO^^cm^^) which belongs to the ISM 
of satellite galaxies is excluded from the measurement. 

3.2 CII Absorption 

In order to more carefully investigate the possibility of de- 
tecting cold filaments using metal absorption lines, we com- 
pute the strength of the CII A1334 absorption. Our choice 
of line is dictated not only by the temperature of the fila- 
mentary gas which is not high enough to significantly pro- 
duce more highly ionized metallic elements such as CIV, but 
also because it is empirically known to yield the strongest 
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Figure 2. An example of the contribution from filamentary gas 
to the CII A1334 absorption. Upper panels show the hydrogen 
column density and mass- weighted metallicity distributions for a 
IO^^Mq halo at 2 ~ 3.8. A cold filament is indicated on the left 
panel, which is receding from the observer and therefore should 
produce a redshifted absorption line. A white square denotes the 
central ~10 kpc^ region over which the absorption spectrum is 
obtained. Bottom panel shows an absorption profile of this galaxy 
(thick grey solid line). We also compute the absorption without 
the central disc by neglecting the opacity of central cells (r < 
0.1r20o) (green solid line) and absorption without the outer gas 
(r > 0.1r20o) (blue dashed line). Also included is the absorption 
produced when the outer gas (r > 0.1r2oo) is assumed to have 
solar metallicity in the absence of the central disc (red dotted 
line). It can be seen that the contribution from the filament to 
the low-ionisation metal line is negligible compared to that of 
the ISM of the galaxy in every case. An uncorrected spectrum 
is included to show the effect of the correction applied to the 
velocity distribution of gas (thin solid grey line). 

absorption feature (jSteidel et al.|[201ol ) We do not attempt 
to model absorption by CII accurately which would re- 
quire detailed radiative transfer, but instead derive an up- 
per limit by making several extreme assumptions. Firstly, 
we assume that all the carbon present in filaments is eligi- 
ble for the CII A1334 transition. Sec ondly, we use the so - 
lar abundance ratio {[C/Z]q ~ 0.178. lAsplund et all [20091 ) 
to obtain the carbon column density for a given metallic- 
ity in the simulation. The optical depth of a grid cell is 
computed as r = crcii^ciiAZ, where ncii is the carbon 
number density, Al is the size of the grid cell, and acii 
is the cross-section for the line transition, which is calcu- 
lated as (Jcii = (37rcrT/8)^/VAo ^ 1-5 x 10"^*cm^ Here 
(Tt is the Thomson cross-section, Ao is the rest-frame wave- 
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length of the transition, and / is the corresponding osciUa- 
tor strength. We then correct the optical depth by assum- 
ing that each grid cell has a Gaussian velocity distribution 
with a dispersion ((7id/2), which is obtained from the line- 
of-sight velocity dispersion (crioa) computed using the closest 
27 neighbouring cells. We have tested the validity of the cor- 
rection by using a high resolution in the Nut series (12 pc 
resolution, Devriendt et al. in prep.), and found that this 
procedure yields a accurate approximation of the maximum 
absorption strength and FWHM of the absorption profile 
that one would derive by using a much higher number of 
resolution elements. Finally stars are assumed to dominate 
the UV emission and we use the iMarastonI (|2005l ') spectral 
energy distributions to derive the continuum flux around 
A ~ 1334A, which depends on their mass, age, and metal- 
licity. The emission from each star particle in the galaxy is 
then used to estimate an observed flux, as attenuated by in- 
tervening gas present along the line of sight. Note that only 
the flux emitted by the central 10 kpc^ of a galaxy is used 
to compute the absorption line so as to mimic the obser- 
vational resolution of FWHM~ 0.40" for galaxies at z ^ 2 
l|Steidel et al.ll2010t ). 

Fig. [2] shows the HI column density map, projected 
metallicity distribution, and the corresponding absorption 
profile for a galaxy residing in a A^200 ~ 10^^ M© halo at 
z ~ 3.8. A dense filament is clearly seen not only in the H 
column density map, but also in the metallicity map, with 
metallicity (-^ 10~'^ — W~^ Zq). The filamentary gas is reced- 
ing from the observer; hence if detectable it should produce 
a redshifted absorption line. However, it turns out that the 
absorption signal is dominated by the ISM of the galaxy 
lighting up the filament. When the absorption due to the 
ISM is arbitrarily removed by neglecting the opacity from 
the gas inside the central gas disc (r < 0.1r2oo), the absorp- 
tion feature vanishes (green line) . Even when the gas outside 
r > 0.1r2oo is assumed to have solar metallicity, the absorp- 
tion strength is still much smaller (red dotted line) than the 
absorption produced by the galaxy's ISM (grey and blue 
dashed lines). This strongly suggests that the primary rea- 
son why it is so difflcult to detect the cold fllament is that 
the density of the fllamentary gas is much lower than that of 
the galaxy's own ISM. As a result, and since the two absorp- 
tion lines are not well enough separated in velocity space, 
the filament absorption signal is completely swamped by the 
high level of ISM absorption in the red wing of the line. 

In order to see if we can bring out the fllamentary signal 
by stacking absorption proflles, we analyse the absorption 
spectra of 132 and 386 massive galaxies (Af2oo > W^^Mq) 
at z — 3.8 and z = 2.5 respectively along 6 projections 
{+x, —X, +y, —y, +z, —z) and flnd that the optical depth for 
the CII A1334 transition by the fllaments is fundamentally 
small regardless of redshift. Fig. [3] shows that the stacked 
(mean) absorption strength is unaffected by the presence of 
fllaments. To check whether the signal from filaments could 
potentially become noticeable if the ISM was less metal- 
enriched, we also examined the case where only a tiny frac- 
tion (1%) of the carbon in the ISM is eligible for the CII 
transition (blue lines in Fig. [3]), but the difference between 
absorption profiles resulting from all the gas along the line of 
sight versus the case where we exclude the gas in the central 
region (r < 0.1r2oo) is still minute. Therefore, we conclude 
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Figure 3. Stacked absorption profiles of massive galaxies 
(M2()o > IQI^Mq) at z = 3.8 (left) and z = 2.5 (right). We 
take 6 projections {+x, —x, +y, —y, +z, —z) for 132 galaxies at 
z = 3.8 and 386 galaxies at z = 1.5. The stacked spectra are 
obtained by taking the mean of the normalised intensity. We also 
show the absorption spectra expected when only a fraction (1%) 
of the ISM carbon is assumed to produce the absorption (blue 
lines, see the text). The contribution from the outer regions (i.e. 
filament if any) is negligible in the absorption profile. 



that the presence of cold filaments is very difficult to confirm 
with low-ionisation metal absorption lines. 

The optical depth of filaments may be under-estimated 
due to the finite resolution of the Horizon-MareNostrum 
simulation. For example, the density of the filamentary 
structure at z = 7 in the Horizon-MareNostrum simu- 
lation (0.005 ^ riH ^ 0.1) is more than an order of magni- 
tude lower than that of the Nut simulation (0.1 mh ;$ 1) 
where the filaments are fully resolved (Powell et al. in prep.). 
However, this increase does not suffice to produce a strong 
absorption signal. Moreover, we find that the effect of in- 
creasing resolution affects the ISM density more consider- 
ably (it becomes more than four orders of magnitude higher 
in the Nut than in the Horizon-MareNostrum simula- 
tion). As a consequence the ISM causes a stronger absorp- 
tion signal, which more than compensates the increased con- 
tribution from the cold filament. 

It should also be noted that finite resolution affects 
the metallicity of the filamentary gas, so that in the real 
Universe its metallicity may be higher than the values (~ 
10"* -10"^^©) we report for the Horizon-MareNostrum 
simulation. Indeed, our simulation cannot resolve all the 
small galaxies that might pollute the pristine filamentary 
gas. Besides, it is well known that the energy from super- 
novae will artificially be dissipated when simulations are run 
with an insufficient level of resolution. Indeed, in such a case, 
supernovae only explode in dense grid cells, resulting in sub- 
stantial radiative losses and negligible momentum transfer 
to the surrounding gas. Under these circumstances, metals 
cannot disperse properly. For comparison, the ultra-high res- 
olution Nut simulation indicates that the metallicity of the 
filamentary gas around a ~ 10^ Mq galaxy can reach val- 
ues up to W~^Zq already at z ~ 7. Yet, such an increase 
in metallicity would still be inconsequential for the absorp- 
tion spectra. Furthermore, we believe that the metallicity in 
the filaments is not likely to rise appreciably beyond these 
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values, because most of the supernova ejecta escapes in a 
direction perpendicular to that of the elongated filament, 
which makes it difficult to efficiently enrich the filamentary 
gas with metals (Geen et al 2010, in prep.). However, for 
the sake of completeness, we present in Fig. [2] the case of 
a halo for which the gas outside r > 0.1r2oo is arbitrarily 
assumed to have solar metallicity. We find that the absorp- 
tion strength of the filament (red dotted line) is still much 
smaller than the absorption produced by the galaxy's ISM. 



4 CONCLUSIONS AND DISCUSSION 

Cosmological simulations predict th at high-z galaxies grow 
by acquiring gas from cold streams (jPekel et al.ll2009l ). but 
no observational confirmation has been obtained yet. Based 
on a statistical sample from the Horizon-MareNostrum 
simulation, we argue that low-ionisation metal absorption 
features, such as CII A1334, arising from intervening cold 
filaments are extremely hard to distinguish from absorption 
by the ISM of high-z star-forming galaxies. This is primarily 
because the optical depth for the low-ionisation transition 
from cold filamentary gas is minuscule, compared with that 
of the ISM of the host galaxy. Moreover, the filamentary 
absorption is not redshifted enough with respect to the ISM 
absorption, so that the residual ISM absorption in the red 
wing of the line is still prominent. This small optical depth 
of filaments mainly finds its source in the intrinsically low 
densities and metallicities of the cold gas when compared to 
those of the galaxy's ISM. Another factor is the geometry of 
the flow which lowers the probability of detecting filaments 
as their column density will rarely be maximised by being 
aligned with the line of sight. 

As an alternative to using a single galaxy, one could 
probe circumgalactic regions by using a paired background 
galaxy. This method alleviates the importance of ISM ab- 
sorption since when probing circumgalactic regions in this 
way, the absorption by the ISM of the foreground galaxy 
will occur at a different spatial position from that pro- 
duced by filaments. Unfortunately, the rare occurrence of 
suitable foreground-background galaxy pairs makes it dif- 
ficult to probe more than one line of sight per foreground 
galaxy. As a result, high resolution individual spectra are 
hard to obtain and o ne has to resort to stacking the spectra 
of multiple galaxies (jSteidel et aLllioiol ). Contrary to what 
might be expected, stacking will wash out the cold filament 
absorption signal since absorption by infiowing gas does not 
neatly separate from that caused by outfiows as was the case 
when probing the circumgalactic region using single galax- 
ies. Indeed, cold filament absorption against the background 
galaxy light will not only be redshifted, but also blueshifted 
as one expects that on average as many cold streams will 
be detected moving towards as away from the observer. Ab- 
sorption by outfiows will also suffer the same fate. This will 
make it all the more difficult to argue whether the observed 
absorption features are driven by infalling gas or by outfiows 
from high-z star-forming galaxies. Moreoever, the metal col- 
umn density of cold filaments is minuscule, as already men- 
tioned. The hydrogen column densities of the cold filaments 
are distributed around lO^^cm"^ at these redshifts, yielding 
corresponding carbon column densities around lO^^cm"^ if 
an average Z = 0.001^© is used. Even if all carbon atoms 



are assumed to be eligible for the CII transition, the optical 
depth is only around 10~^ in the line. On the other hand, 
outflows are expected to be very metal rich, so that even 
though higher transitions like CIV are expected to domi- 
nate the absorption signal, the amount of CII absorption 
from these outflows might still swamp that produced by the 
cold filaments. 

Finally, we have assessed possible numerical resolution 
issues on column density and metallicity of the filaments 
using the very high resolution numerical simulation Nut 
suite and found that our conclusions remain by and large 
unchanged. 

Based on these considerations, we conclude that the 
presence of the cold filament is difficult to disprove/prove 
with low-ionisation metal line absorption. Instead, the Ly- 
man Q emission route seems more promising to detect cold 
filaments, but the line profile is more sensitive to the kine - 
matics of the intervening gas (|Verhamme et ai]|2006l . [2008l V 
Therefore this will require full blown radia tive transfer cal- 
culations (e.g. iFaucher-Giguere et al.ll20ld l. 
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